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SECTION  1 
INTRODUCTION 


This  report  documents  modifications  and  additions  incorporated  into  the 
ABRES  Shape  Change  Code  (ASCC),  originally  u?velooed  under  the  Passive  Nosetip 
Technology  (PANT  II)  program  (Reference  1),  and  subsequently  modified  and 
improved  as  part  of  the  Reentry  Vehicle  Technology  (REV-TECH)  program 
(References  2  and  3).  In  this  document,  ASCC 7 7  denotes  the  original  version 
of  ASCC  developed  under  PANT  II,  and  the  current  version  of  ASCC  is  referred 
to  as  ASCC80. 

The  overall  objectives  of  the  Aerodynamic  Heating  Computations  for 
Projectiles  program  were  threefold: 

1.  Modify  the  in-depth  heat  conduction  package  to  improve  ASCC's 
capabilities  to  handle  slender  multimaterial  configurations 

2.  Extend  the  developments  of  planar  ASCC  modifications  to  predict 
heating  of  swept  fin  configurations  to  include:  (a)  turbulent  flow 
on  swept  wings;  (b)  2-D  shock  shape;  and  (c)  improved  in-depth  heat 
conduction  routines 

3.  Develop  an  interactive  computational  grid  developing  routine  to 
simplify  the  procedure  for  inputting  body  configurations  and 
developing  computational  grids 

The  modifications  made  to  ASCC80  covering  the  first  objective  are 
documented  in  Volume  I  of  this  report.  Changes  made  to  ASCC80  covering 
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Objectives  2  and  3  are  documented  in  Volumes  II  and  III,  respectively.  In 
this  document,  the  updated  ASC  code  is  referred  to  as  BRLASCC. 

8RLASCC  includes  all  the  modifications  and  improved  capabilities  made 
to  ASCC  under  the  previously  mentioned  programs  and  the  following  additional 
modi fications: 

•  Body  points  which  are  positioned  behind  the  origin  of  rays  now  have 
their  rays  perpendicular  to  the  axis  of  symmetry,  rather  than 
emanating  from  the  single  origin.  This  extends  the  code 
applicability  and  accuracy  for  more  slender  nosetips,  precluding 
excessive  "skewing"  of  the  mesh. 

•  The  user  can  now  vary  the  implicit  grid  thickness  from  ray  to  ray 

•  The  interface  modeling  between  different  materials  has  been 
improved,  particularly  in  the  explicit  grid,  so  that  material 
properties  used  in  the  finite-difference  conduction  equations  are 
utilized  more  accurately 

■»  The  user  may  now  input  latent  heat  of  fusion  for  melting  materials 
and  BRLASCC  will  account  for  this  energy  during  melting  in  the 
in-depth  conduction  solution 

®  The  thermal  resistance  in  the  finite-difference  conduction 

equations  for  the  implicit  and  explicit  grid  have  been  modified  to 
include  contact  resistances,  which  will  be  input  by  the  user 

•  Up  to  six  materials  may  now  be  input  with  up  to  five  allowable 
interfaces  along  each  ray  within  the  implicit  grid 

•  The  common  blocks  nave  been  restructured  to  allow  individual  users 
to  easily  identify  and  specify  implicit  and  explicit  grid 
parameters  and  recompile  the  code  when  necessary  using  the  VAX 
Fortran  "INCLUDE"  statement 
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Technical  discussion  of  the  BRLASCC  modifications,  including  rationale 
is  presented  in  Section  2.  Section  3  is  devoted  to  a  discussion  of  input  and 
output.  For  user  convenience,  a  complete  set  of  input  instructions  is 
included.  The  user  is  encouraged  to  refer  to  Reference  2  for  contrast  of  the 
input  requirements. 
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SECTION  2 

TECHNICAL  DISCUSSION 


Significant  modifications  and  improvements  have  been  incorporated  in 
ASCC80  to  generate  BRLASCC.  As  stated  in  the  introduction,  the  changes  we^e 
aimed  largely  at  improving  ASCC's  capabilities  to  handle  the  longer,  more 
slender  multimaterial  configurations  of  interest  to  the  U.S.  Army  Ballistic 
Research  Laboratory.  BRLASCC  includes  the  following  improvements: 

(1)  accommodation  of  the  entire  projectile  configuration  by  removing 
difficulties  which  result  from  extreme  skewness  of  the  implicit  rays  with  the 
surface,  (2)  improved  thermal  resistance  modeling,  including  the  inclusion  of 
contact  resistance  to  the  finite-difference  conduction  equations,  (3)  improved 
in-depth  modeling  by  inclusion  of  latent  heat  of  fusion,  (4)  increased 
capability  to  handle  up  to  six  different  materials  and  allowing  up  to  five 
interfaces  on  a  ray  within  the  implicit  layer,  (5)  restructured  common  blocks, 
allowing  the  user  flexibility  for  modifying  the  array  sizes  which  determine 
the  limits  for  the  implicit  and  explicit  grids.  A  discussion  of  these 
modifications  is  presented  in  Sections  2.1  through  2.5. 

2.1  ACCOMMODATION  OF  THE  ENTIRE  PROJECTILE  CONFIGURATION 

ASCC80  utilizes  dual,  overlapping,  orthogonal  coordinate  systems 
composed  of  the  following:  (a)  moving  body-oriented  coordinate  system, 
(s,r,v),  which  extends  over  the  heated  surface  layer,  and  (b)  a  fixed 
cylindrical  coordinate  system,  (x,y,y),  extending  over  the  entire 
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computational  domain  (Figure  2-1).  Different  finite  difference  equations  are 
used  in  each  region:  implicit  in  (a)  and  explicit  in  (b).  Therefore,  the 
moving  coordinate  system  is  referred  to  as  the  implicit  grid  or  surface  layer, 
and  the  fixed  coordinate  system  is  referred  to  as  the  explicit  grid.  The 
major  difficulty  associated  with  the  use  of  this  gnddiny  scheme  is  the  use  of 
a  single  origin  of  rays  for  the  implicit  grid.  This  procedure  was  originally 
developed  and  validated  for  nosetip  problems,  however  two  potential 
computational  problems  arise  in  the  conduction  solution  when  used  for  .  ie 
prediction  of  more  slender  vehicle  configurations  such  as  nosetip  and 
heatshield. 

The  first  problem  pertains  to  the  placement  of  the  origin  of  rays,  OX. 

A  basic  assumption  in  the  formulation  of  the  finite-difference  equations  for 
the  implicit  layer  is  that  the  rays  be  as  nearly  ortho- normal  to  the  surface 
as  possible.  ASCC  contains  approximate  corrections  to  account  for  small 
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Figure  2-1.  In-Depth  Conduction  Coordinate  System  and 
Finite-Difference  Grid 
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nonorthogonalities,  however  these  corrections  are  inadequate  for  cases  where 
deviation  from  the  assumption  of  orthogonality  is  significant.  On  long, 
slender  configurations  it  is  Impossible  to  meet  this  criterion.  If  OX  is 
placed  near  the  nosetip,  the  implicit  layer  becomes  so  thin  and  distorted  on 
the  heatshield  that  the  assumption  of  near  orthogonality  is  invalid.  If  OX  is 
placed  far  back  on  the  frustum,  the  surface  layer  near  the  tangent  point  has 
the  potential  of  forming  a  cusp  in  the  implicit  grid,  causing  instabilities  in 
the  conduction  solution.  This  results  in  negative  temperatures  and  negative 
time  steps  which  halt  further  execution  of  the  code.  By  judicious  placement 
of  the  origin  OX  and  using  proper  mesh  size  for  the  explicit  grid,  reliable 
predictions  can  be  made  on  the  entire  configuration.  However,  the 
computational  speed  is  slow  due  to  small  time  steps  required  by  the  explicit 
procedure  stability  criterion,  which  is  affected  by  the  smaller  spatial  mesh 
size. 

The  second  problem  arises  due  to  the  limitations  placed  on  the  input  of 
implicit  layer  thickness  along  the  rays.  This  problem  is  illustrated  in 
Figure  2-2  and  i$  caused  by  the  constant  thickness  of  the  implicit  layer  along 
each  ray.  Surface  Layer  1  of  Figure  2-2  uses  a  thic1,  layer  which  strongly 
skews  the  back  of  the  implicit  mesh  near  the  origin,  resulting  i,i  the  cusp. 
This  distortion  of  the  mesh  can  cause  numerical  instabilities  resulting  in 
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Figure-  2-2.  Skewing  of  Implicit  Grid  Resulting  From 
Thick  Surface  Layer 
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negative  temperatures.  Surface  Layer  2  uses  a  thinner  layer,  and  the 
distortion  is  nearly  eliminated  in  the  vicinity  of  the  origin;  however, the 
layer's  thickness  with  respect  to  the  explicit  grid  y-direction  is  very  small 
requiring  small  explicit  mesh  spacing  and  resulting  in  small  time  step  sizes. 

The  limitations  on  the  location  of  the  origin  of  rays  and  the  implicit 
layer  thickness  poses  restrictions  on  the  mesh  si/e  of  the  explicit  grid,  ing 
consequently  the  length  of  the  vehicle  which  can  be  accurately  modeled.  Tg 
alleviate  these  difficulties  in  BRLASCC,  two  modifications  were  made. 

First,  the  old  “radial"  implicit  gridding  scheme  extant  in  ASCC80  has  r.een 
converted  to  a  dual  system  which  is  “radial"  in  the  nosetip  region  an'; 
"cylindrical"  in  the  afterbody  region.  The  difference  between  the  old  and  r.-- 
grid  systems  is  illustrated  in  Figures  2-3  and  2-4.  For  body  points  which  1; 
behind  the  origin  of  rays,  the  rays  are  new  perpendicular  to  the  axis  of 
symmetry  instead  of  originating  at  OX.  This  will  allow  the  user  to  use  the 
ray  system  for  shape  change  with  in-depth  conduction  in  the  nose  region  and 
will  prevent  the  severe  nonorthogoral i ty  with  tne  surface  at  distances 
downstream  of  the  nosetip. 

Secondly,  the  current  requirement  that  the  implicit  grid  thickness  be 
constant  along  each  ray  has  been  removed.  Instead  of  inputting  a  table  of 
implicit  grid  node  spacing,  the  code  now  accepts  a  table  of  normalized  node 
spacing  and  a  table  of  implicit  layer  total  thickness  along  each  ray.  This 
will  allow  the  user  to  use  thicker  and  better  contoured  implicit  layers  as 
seen  in  Figure  2-4  and  avoid  the  skewing  shown  in  Figure  2-2.  On  the 
frustum,  many  more  explicit  nodes  fall  within  the  implicit  layer  giving  tette 
definition  to  the  explicit  boundary  conditions.  The  skewing  near  OX  is 
eliminated  through  proper  use  of  the  implicit  grid  thickness  table. 

Figure  2-5  compares  in-depth  temperature  prediction;  made  using  ALCC8U  and 
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(b)  Conduct  ton  Grid  tHIno  th«  Wod'Mdd  ASCC  SctidM 


-4.  Conduction  Grid  Using  BRLASCC  Scheme 
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Figure  2-5.  Comparison  of  In-Depth  Temperature  Profiles  Predicted  With 
ASCC80  and  BRLASCC  Gridding  Schemes 

BRLASCC  gridding  schemes.  As  can  be  seen,  BRLASCC  provides  much  greater 
coupling  between  the  explicit  (in-depth)  temperatures  and  the  implicit 
(surface  layer)  temperatures. 

2.2  IMPROVED  THERMAL  RESISTANCE  MODELING 

The  conduction  package  in  ASCC  allows  several  material  interfaces  to  be 
included  in  the  analysis.  Each  expli.it  grid  point  (node)  and  implicit  grid 
point  (nodlet)  is  flagged  with  an  inteqer  denoting  the  material  number  of  that 
grid  point.  ASCC80  uses  a  simple  average  of  the  material  property  values  of 
the  node  materials  to  model  the  properties  of  each  finite-difference  cell  and 
does  not  account  for  the  contact  resistance  at  material  interfaces.  When 
adjacent  materials  have  very  dissimilar  physical  properties  this  assumption 
can  cause  substantial  smoothing  and  loss  of  accuracy.  In  addition,  ASCC80 
contains  no  provision  for  the  inclusion  of  contact  resistance  between 
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materials.  BRLASCC  overcomes  these  problems  by  calculating  the  proper  thermal 
resistance  between  nodes  and  keeps  track  of  interface  locations  In  both  the 
implicit  and  explicit  grids  and  Includes  contact  resistance  If  input  by  the 
user. 

2.2.1  Implicit  Grid  Modifications 

The  conduction  equation  in  the  moving  orthogonal  coordinate  system 
(implicit  grid)  under  the  axisymmetric  assumption  (a/3r  *  0)  is: 


r  3T 
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where 

Cp  =  specific  heat 

r0  =  body  circumferential  radius  of  curvature 
rb  *  ro  +  r*cos(e) 

rc  =  local  streamwise  radius  of  curvature 
•c  =  thermal  conductivity 
p  =  density 

n  =  surface  normal  recession  rate,  n  *  -r 
T  =  temperature 
t  =  time 

9  =  angle  between  normal  to  local  surface  and  axis  of  symmetry 
s  =  streamwise  distance  along  body 

r  =  distance  normal  to  body  surface  at  s,  measured  from  the  surface 


16 


where  is  the  body  radius  of  curvature  in  a  plane  perpendicular  to  the  axis 
rc  is  the  local  surface  streamwise  radius  of  curvature,  and  e  is  the  surface 
inclination  with  respect  to  the  axis. 

The  finite-difference  equation  for  the  implicit  grid  in  ASCC80  is: 
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where  the  superscript  prime  indicates  the  quantities  evaluated  at  the 
new  time,  q-|c  is  the  lateral  conduction  correction  term  associated  with 
nonorthogonal  rays,  and  qs  is  the  lateral  diffusion  term.  In  the  formulation 
of  the  finite-difference  expression,  it  is  assumed  that  the  lateral  diffusion 
term  is  small  compared  with  the  normal  diffusion  term  and  therefore  may  be 
evaluated  in  terms  of  the  old  temperatures  and  appears  as  a  "source"  term  in 
the  equation. 
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The  finite-difference  equation  for  the  implicit  grid  has  been  modified 


and  appears  in  BRLASCC  as: 
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where  the  thermal  resistances  are  given  as: 
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Note  that  Instead  of  the  simple  average  of  thermal  conductivity  used 
in  ASCC80,  8RLASCC  calculates  the  thermal  resistance  through  each  material  as 
a  function  of  conduction  path  length  and  thermal  conductivity,  and  includes 
contact  resistance. 

2.2.2  Explicit  Grid  Modifications 

The  conduction  equation  in  the  fixed  cylindrical  coordinate  system 
(explicit  grid)  is: 
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The  finite-difference  equation  for  the  explicit  grid  in  ASCC80  is: 
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1+1/2 ,0 


(TI,0-  T1 


(3) 


where  again  the  superscript  prime  indicates  the  quantity  to  be  evaluated  at 
the  new  time.  The  thermal  resistance  is  given  as: 


(6) 
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The  thermal  resistance  term  has  been  modified  In  BRLASCC  and  appears 


as: 


R 


I.J-1/2 


8y  '  »■»-!  .  '  V  i 


\)-l 


contact 


;  6  i  distance  from  y,  j  (7) 
*  to  Interface  inJ~ 
y-di rection 


lI  +  l/2,J 


Vl  -  5x  .  5x  -  XI  .  1 

~ ' - * — ’ — 


kI+l 


'I 


contact 


;  6  =  distance  from  x^ 

x  to  interface  in1 
x-di rection 


(8) 


As  in  the  expression  for  the  implicit  conduction  equation,  the 
thermal  resistance  calculation  has  been  modified  to  use  the  appropriate 
materia1,  thermal  conductivity  and  conduction  path  length  and  includes 
contact  resistance. 

2.3  LATENT  HEAT  OF  FUSION 

A  model  to  account  for  the  latent  heat  of  fusion  during  in-depth 
melting  of  a  material  has  been  included  in  8RLASCC.  The  energy  storage  terms, 
pCp,  are  affected  in  both  the  implicit  and  explicit  formulations  of  the 
conduction  equation.  The  energy  required  to  melt  a  material  may  now  be  input 
to  the  code  by  the  user  (see  Section  3,  Input  Table  6  —  Material  Properties) 
and  is  included  in  the  sensible  enthalpy  calculation  for  the  material  tables. 

A  saw-tooth  function  of  specific  heat  versus  temperature  Is  constructed  from 
user  input  of  latent  heat  of  fusion  (A L),  melt  temperature,  specific  heat  of 
the  solid  at  melt  temperature,  specific  heat  of  the  liquid  at  melt 
temperature,  and  the  temperature  difference  over  which  the  user  would  like 
melt  to  occur  (ATm).  The  integration  of  Cp  over  AT  results  In  the  proper 
sensible  enthalpy,  as  shown  in  Figure  2-6. 


2G 


Figure  2-6.  Saw-Tooth  Function  of  Specific  Heat  Versus  Temperature 
for  a  Melting  Material 


An  analytical  solution  is  available  from  Carslaw  and  Jaeger 
(Reference  4)  for  the  melting  of  a  semi-infinite  solid  without  melt  removal. 

In  this  solution  the  transient  temperature  distribution  and  position  of  melt 
front  are  given  as  a  function  of  time,  material  properties,  and  initial 
condition.  This  solution  was  used  as  a  reference  to  check  against  the  8RLASCC 
solution.  Three  cases  were  run: 

•  Implicit  layer  thickness  of  1  inch;  uniform  implicit  nodlet 
spacing  of  0.14286  inch;  uniform  explict  node  spacing  of 
0.50  inch 

•  Implicit  layer  thickness  of  1  inch;  uniform  implicit  nodlet 
spacing  of  0.07143  inch  (half  of  case  1);  uniform  explicit  node 
spacing  of  0.50  inch 
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•  Implicit  layer  thickness  of  1  inch;  uniform  implicit  nodlet  spacing 
of  0.07143  inch;  uniform  explicit  node  spacing  of  0.25  inch  (half 
of  case  1) 

The  material  properties  used  in  the  analysis  are  given  in  Table  2-1. 

The  Carslaw  and  Jaeger  solution  and  BRLASCC  solutions  are  compared  in 
Figures  2-7  and  2-8  for  temperature  distribution  at  200  seconds  and  melt 
front  position  versus  time,  respectively.  The  results  show  good  agreement 
between  BRLASCC  and  the  analytical  solution.  However,  it  is  clear  from  the 
figures  that  the  accuracy  of  the  solution  is  mesh  si2e  dependent;  halving 
the  mesh  spacing  nearly  halves  the  error.  The  user  must  exercise  good 
judgment  in  mesh  spacing  when  generating  solutions  involving  in-depth 
melting. 

The  user  must  also  note  that  the  model  for  latent  heat  of  fusion  is 
not  adequate  for  surface  melting  since  melt  removal  coupled  with  shape 
change  is  not  modeled.  The  surface  energy  balance  would  be  in  error  as  melt 
occurs,  resulting  in  incorrect  shaping,  and  these  errors  would  grow  in  time. 

Table  2-1.  Melting  of  a  Semi -Inf inite  Solid  —  Material  Properties 


Material  Properties 

Solid 

Liquid 

k  (Btu/sec-ft-°R) 
a  ( ft'/sec) 

Cp  (Btu/lbm-°R) 

At  =  117  Btu/lbm 

T0  =  500°F 

Tj  =  2 ,800°F 

V  =  4,000°F 

7.56  x  lO-3 

1.208  x  10-4 
0.128 

Latent  hea 
Initial  soli 
Melt  te 
Surface  t 

2.10  x  10-3 
2.222  x  10-5 
0.1932 
t  of  fusion 
d  temperature 
mperature 
emperature 
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Tetnperatu 


Distance  Measured  from  Surface,  X  (In) 


Ficrjre  2-7.  Comparison  of  In-Depth  Temperature  Distribution 

for  a  Melting  Solid  Predicted  by  Carslaw  and  Jaeger 
and  BRLASCC 


Melting  of  a  Semi  -  Inf ini te  Solid 


□  □ 


G 


■  \  Unitor*  mplicil  fcrid 
(pacing  •  0  07M3  m 
Unitor*  eiplicii  grid 
(pacing  ■  0  2$  in 

_J _ 1 - J - 1 - 

50  100  150  .  200  2! 

Time  (sec) 


2-8.  Comparison  of  Melt  Front  Location  for  a  Melting 
Solid  Predicted  by  Cars! aw  and  Jaeger  and  BRLASCC 


2.4  MULTIPLE  MATERIAL  HANDLING 


As  the  application  of  the  ASC  code  broadens,  more  complex  internal 
geometries  are  being  considered.  ASCC80  contains  an  easily  used  general 
shape  input  option  that  allows  the  description  of  complex  internal 
geometries  with  up  to  four  materials,  and  allowing  two  interfaces  on  each 
ray  within  the  implicit  layer.  BRLA5CC  has  been  updated  to  allow  the  input  of 
up  to  six  materials  and  allows  up  to  five  interfaces  on  each  ray  within  the 
implicit  layer.  This  will  facilitate  the  calculation  of  complex  in-depth 
configurations. 

2.5  FLEXIBLE  GRID  ARRAY  SIZES 

The  size  of  the  arrays  used  for  the  thermal  conduction  grids  is  often 
restrictive  for  complete  configurations  of  interest.  This  is  especially 
true  of  the  number  of  y-values  in  the  explicit  grid.  Increasing  the  alloweo 
grid  sizes  is  simply  a  matter  of  increasing  the  dimensions  of  the  arrays. 

Since  BRLASCC  makes  extensive  use  of  the  VAX  Fortran  INCLUDE  statement,  this 
change  only  requires  modifications  to  the  common  blocks  containing  the 
necessary  arrays  and  a  recompilation  of  the  code.  It  should  be  noted  that 
redimensioning  of  those  arrays  used  in  the  restart  routine  DATFIL  would  be 
necessary  to  preserve  the  restart  capability. 

One  difficulty  with  increasing  the  dimensions  of  the  arrays  will  be 
the  core  limitat’ons  of  the  oarticular  machine  being  used.  For  this  reason, 
the  arrays  to  be  redimensioned  have  been  regrouped  in  their  common  blocks 
and  INCLUDE  files  to  facilitate  changing  the  grid  array  sizes  whenever 
necessary.  These  groupings  of  common  blocks  are  commented  for  user 
convenience.  The  array  sizes  have  not  been  changed  between  ASCC80  and 
BRLASCC;  this  is  left  to  the  discretion  of  the  user.  A  sample  of  the 
regrouped  and  commented  common  blocks  for  use  with  the  INCLUDE  statement  is 


shown  in  Table  2-2. 


ooooooo  o  oooooo 


Table  2-2 


A  Sample  of  Commented  Common  Blocks 


*********************************************************************** 

THE  FOLLOWING  COMMON  BLOCKS  CONTAIN  CONSTANT  VARIABLES  WHOSE 
VALUES  WILL  NOT  CHANGE  DURING  EXECUTION 
*********************************************************************** 


COMMON  /CONST/ 


*  ALMAX, 

ALMIN, 

OEXMIN, 

OLTMIN, 

DLTRAN, 

*  DOX, 

DSHANG, 

FMO, 

GFHI, 

HJNCTN, 

*  HSCF, 

HSCI , 

HTJF , 

HTJI , 

HZEROT, 

*  OX, 

OY, 

PR, 

*  RBAS, 

RNI , 

RPHI , 

RSIDE , 

RX, 

*  RZ. 

SPHTC, 

STRD , 

THETA, 

THETAC, 

*  TRNTIM, 

VLN , 

WT, 

XCG, 

ZMAX, 

*  ZMAXTJ , 

ZSIDE1 , 

ZSTAG1 , 

*  KM ( 30 ) , 

REM( 30) , 

*  REC0SD(36) 

*  PRTBL (60) , 

PRTBl{60) , 

*  TIMUSR(?50) 

COMMON  /CONST I / 

*  I ATM, 

I ATMS, 

IBRUPT, 

ICARB, 

iorctn. 

*  I FACE, 

IFLG09 , 

IHMAX, 

IE , 

INDATM, 

*  IMOD, 

INPUT, 

IPMAX , 

*  IPRFLG, 

IPRNT , 

IRON, 

IRSTRT, 

ISHFLG, 

*  ISS, 

JL, 

LG, 

EL, 

MATN , 

*  MLAYRS, 

NAM, 

NCL , 

*  NIF, 

NLAYRS, 

NOHEAL , 

*  NOSLO, 

NPRTBL , 

NREYCR , 

NSHTBL, 

NTFIX, 

*  NTIMT , 

NTMUSR 

DIMENSIONED 

FOR  IS  IMPLICIT  NODLETS 

,  50  SURFACE 

POINTS, 

X ( 60 )  BY  Y(25)  EXPLICIT 

NODES,  AND 

100  ENTRIES 

IN  GENERAL 

INTERFACE 

TABLES 

COMMON  /CONGRD/ 


*  DSMOV , 

OTOROP,  OXN, 

*  DELN( 15 ) , 

*  Y(25), 

YD I F ( 25 ) , 

*  OEL (50) , 

XI  NIT (50 ) ,  Y INIT(50) 

*  X  ( 60 ) , 

XDIF (60) , 

*  ZIS(IOO), 

RIS(IOO) , 

*  DMIN(60 ,25) 

• 

*  I DROP, 

IMOVE, 

*  NBS(IOO), 

*  NMAT(6G,25) 
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Table  2-2.  A  Sample  of  Commented  Common  Blocks  (Concluded) 


****★★■*★**  +  +  +  Hr  ***********************  <r  ******  ***  *★#*★*  ************* 

THE  FOLLOWING  COMMON  BLOCKS  CONTAIN  VARIABLES  WHOSE  VALUES 
WILL  CHANGE  DURING  THE  EXECUTION 

** ********************************************************************* 


COMMON  /ENVR/ 


*  ALTINF, 

AMACH, 

Al, 

A2, 

CDNT , 

*  DEN, 

DUDZ, 

EMWT2 , 

EMW2 , 

EMW1 , 

*  E2, 

GAM1 , 

GAM2, 

HETAUG, 

HT2 , 

*  HI, 

H2, 

PT2 , 

PI, 

P2, 

*  REYCR, 

RN, 

R0T2 , 

R01 , 

R02, 

*  SSONIC , 

STRAN, 

TSTAGP , 

TT2, 

Tl, 

*  T2, 

UR  1 , 

VIST2 , 

V I S2 , 

VI , 

*  V2, 

*  IESTAT, 

IEXX  , 

INOSE, 

KSHOLD , 

LCT , 

*  ltt. 

NOSTRN 

,  NPGENV, 

NT, 

NTS, 

*  NTT 

DIMENSIONED 

FOR  50 

SURFACE  POINTS, 

5  MATERIAL 

INTERFACES 

AND  15  IMPLICIT  NODLETS 


COMMON  /RECSA/ 

*  BPSP(50) ,  CMDX(50)  , 

*  DFIF(50) ,  DPART(50)  , 

*  FVW( 50)  ,  GKR(50), 

*  RSP(50)  ,  RSPNUI 50)  , 

*  SP(50),  SR AY ( 50 )  , 

*  VIMP(50)  ,  ZSP(50)  , 

*  BLEN(50 ,5)  ,  RI(50,5)  , 

*  IHI (38)  ,  ILO( 38)  , 

*  NB( 50 ,6)  , 

*  IMAT(50,15) 


CMFX( 50)  ,  DELKE(50), 

EFFK(50)  ,  EMD0T(50),  FH50)  , 

HRSP( 50)  ,  PRESP( 50)  , 

RUCHSP(50),  SDOT ( SO )  ,  SD0TE(50)  , 

TANF  1(50)  ,  THETSP( 50)  ,  TSP(50), 

ZSPNU( 50)  , 

ZI(50 ,5) , 

IR( 38)  , 
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SECTION  3 


INPUT  AND  OUTPUT 


The  following  section  is  a  user's  guide  of  the  input  requirements  for 
the  modified  ABP.ES  Shape  Change  Cod?  (BRLASCC).  The  input  has  changed  from 
ASCC80  in  the  following  manner: 

•  On  the  first  restart  card,  variable  OOX  has  been  added.  When  the 
automatic  shifting  of  origin  option  is  exercised,  BRLASCC  will  shift 
the  origin  by  OOX  inches  rather  than  OX  inches  as  in  ASCC80 , 

•  In  Input  Table  3,  the  ASCC77  general  body  shape  input  option  h<.s  been 
eliminated.  All  general  shapes  and  multimaterial  cases  must  be 
input  with  the  ASCC80  general  shape  input  methodology. 

•  In  Input  Table  3,  variable  implicit  layer  thickness  on  each  ray  is 
now  an  option.  The  user  must  specify  the  normalized  spacing  tnat 
will  apoly  to  every  ray  and  the  implicit  layer  thickness  on  each 
ray. 

•  In  Input  Table  6,  material  properties  for  melting  materials  is  now  an 
input  option.  Should  the  user  wish  to  include  in-depth  meltiny  of 
materials  the  user  must  specify  the  latent  heat  of  fusion  (XLATNT), 
melt  temperature  (TMEIT),  temperature  difference  over  which  melt  is 
to  occur  (OTMELT),  specific  heat  of  the  solid  at  TMELT  (CPSOLD),  and 
the  specific  heat  of  the  liquid  at  TMEIT  (CPLIQO). 
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3.1  INPUT  INSTRUCTIONS 


The  Input  to  the  code  can  be  read  from  a  disk  file  or  data  cards  for  an 
Initial  run,  or  from  magnetic  tape  or  disk  for  a  restart  run.  The  details  of 
the  Input  for  each  of  these  types  of  runs  are  described  below.  The  basic 
Input  for  an  initial  run  consists  of: 

•  Two  restart  Information  cards 

•  Three  title  cards 

•  Nine  input  tables 

Not  all  nine  of  the  Input  tables  are  required  for  every  r  « .  Each 
table  is  preceded  by  a  single  card  containing  the  identifying  te  .<*  umber. 

For  a  restart  run,  only  the  restart  information  cards  are  required, 
and  the  rest  of  the  data  is  read  from  magnetic  tape  or  disk. 

Tne  following  are  FORTRAN  I/O  unit  assignments. 


Unit  Device  Type 


Purpose 


1  Terminal  (Interactive) 
Disk  file  (tatch) 


Lists  current  environment 
number,  time,  and  altitude. 


5  Disk  f’le  (interactive/batch)  Input  data 


6  Disk  file  (interactive/batch)  Output  data 

Printer  (batch) 

Terminal  (interactive) 


14  Disk  file  (Interactive/batch)  Restart  file 

Magnetic  tape  (interactive/batch) 


20  Disk  file  (Interactive/batch)  For  coupled  trajectory  option, 

trajectory  scratch  file 


The  following  sections  describe  the  restart  Information,  title  ca.ds, 
and  the  nine  input  tables,  respectively. 

Restart  Information 

The  "restart"  cards  are  the  first  card  set  in  the  data  deck.  If  the 
run  is  a  restart,  they  are  the  only  cards  required.  These  ct'rds  tell  the  code 
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the  environment  number  from  which  to  begin  the  restart  and  what  modifications 
should  be  made  In  the  conduction  grid.  For  an  initial  run,  the  IRSTRT  field 
should  be  left  blank.  The  BRIASCC  code  stores  the  contents  of  its  common 
blocks  on  the  restart  file  at  the  end  of  each  cycle  of  environment  and/or 
shape  change  calculation.  All  restart  reading  and/or  writing  is  done  on 
Logical  Unit  14  and  it  should  be  assigned  accordingly. 

Two  options  allow  modifications  to  the  implicit  conduction  grid  during 
a  calculation.  These  modifications  may  be  done  automatiCclly  or  by  restarting 
the  solution.  The  first  option  drops  surface  points  from  the  grid  and  allows 
larger  time  steps.  The  second  option  shifts  the  origin  of  rays  and 
interpolates  for  temperatures  at  the  new  implicit  grid  points.  The 
modifications  are  done  automatically  by  the  code  when  the  user  specifies  the 
proper  values  of  IDROP,  IMOVE,  DOX,  DTDROP,  DSMOVE.  They  can  be  done  manually 
by  restarting  the  solution  and  specifying  OXN  and  1PDRP. 

The  format  for  the  restart  cards  is  as  follows. 

Card 


No. 

Columns 

Format 

Data 

Units 

1 

1-  5 

15 

IRSTRT  —  ENVIRONMENT  NO.  to  restart  from 

-- 

6-10 

15 

NTFIX  —  ENVIRONMENT  NO.  to  fix  shape. 

This  causes  the  shape  to  remain  fixed 
subsequent  to  assigned  environment  nunber 

— 

★ 

11-15 

15 

IDROP  —  Flag  to  al^ow  automatic  dropping 
of  surface  points 

-- 

s  0  —  Surface  points  will  not  be  dropped 

3  1  —  Surface  points  will  be  dropped  when 
the  lateral  conduction  time  step 
becomes  less  than  DTDROP 
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Card 

No.  Col umns  Format  Data  Units 

*  16-20  15  IMOVE  —  Flag  to  allow  automatic  shifting 

or  origin  of  rays 

=  0  —  Origin  of  rays  will  not  be  shifted 

=1  —  Origin  of  rays  will  be  shifted  DOX 
inches  when  the  implicit  grid  is 
less  than  DSMOV  inches  from  origin 


21-32 

E12.5 

OXN  --  Z-coordinate  for  the  new  origin  of 
rays  during  restart 

Incn 

33-44 

E12.5 

DOX  --  Origin  of  rays  will  be  shifted 
when  needed  by  this  amount;  1M0VE  =  1 

: .n 

45-56 

E12.5 

DSMOV  —  See  IMOVE 

Inch 

57-68 

E12.5 

DTOROP  —  See  IDROP 

Sec 

1-50 

5011 

-  -  -  Regrid  option  only  ----------- 

(Read  only  if  OXN  *  0.0) 

IPDRP(K)  —  For  K  =  1.  KLF  (50  to  a  card), 
flag  for  dropping  surface  point 

=  0  --  Retain  point 
=  1  --  Drop  point 


Title  Information 

The  next  three  cards  contain  title  information  in  Columns  1  through  72. 
Contents  of  Columns  61  through  72  on  Card  3  are  printed  on  every  page  of  the 
output . 
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INPUT  TABLE  1.  GENERAL  PROGRAM  CONSTANTS  AND  TIME  INFORMATION 


This  table  supplies  the  code  with  computation 


time  information  and 


program  flags  which  indicate  options  to  be  subsequently  read. 


Card 

No.  Col umns  Format  Data 


1  1-  2 

2  3-14 

15-26 

27-38 


3  1-  5 


♦ 


6-10 


12  Enter  U1  (table  number) 

E12.5  —  Initial  time  (for  LG  =  -1  it 

corresponds  to  the  initial  conditions 
for  trajectory  calculations) 


E12.5  TF  —  Final  time 

E12.5  DLTM1N  --  Minimum  allowable  time  step.  If 
negative,  only  user's  specified  time  steps 
10LTT)  are  used.  A  default  value  of  10-® 
sec  will  be  used  if  zero  is  entered.  This 
default  value  is  satisfactory  in  most 
cases;  thus,  the  user  can  usually  enter 
zero 


15 


15 


-o 

3 

c r 
<u 

CSI 

a> 

xi 


Environment  flag 

Flight  (internally  computed 
trajectory.  Table  2  not  needed) 

Fl ight 

Wind  tunnel 

Ballistic  range 

General 

Arc  heater 

Shape  change  flag 


0  —  Shape  change  with  transient  in-depth 
conduction 

1  --  Shape  change  with  steady-state 

in-depth  conduction 

2  --  No  shape  change  (boundary  layer 

solution  only) 


Units 


Sec 

Sec 


Sec 


INPUT  TABLE  1.  Continued 


Card 

No. 


Columns  Format 


11-15 


16-20 


15 


IS 


T3 

01 


z> 

cr 


O 

c 


o 


£> 

<15 


Data 

IPRNT  —  Output  print  flag 
0  --  Summary  table  only 

1  --  Detailed  output  at  body  points 

2  —  Detailed  output  at  integration  points 

Positive  IPRNT  —  Output  at  user 
specified  environment  times  only  (DLTT) 

Negative  IPRNT  --  Output  at  all 
environment  times 

NREVCR  —  Transition  criteria  flag 

1  --  Laminar  flow  (Table  04  not  required) 

2  --  Critical  momentum  thickness  Reynolds 

number  versus  edge  Mach  number 

3  Critical  stream  length  Reynolds 
number  versus  edge  Mach  number 

4  --  Nosetip  axial  distance  versus: 

--  Altitude  for  LG  =  1 
--  Time  for  LG  =  2,  3,  4,  5 
b  —  Anderson  nose  criterion 

\0*7  ZSS  for  onset 

'9  iij  215  for  location 
=  Tw/Te 

LORN  cone  criterion 
6  —  Anderson  nose  criterion 

\P*7  255  for  onset 


k  1 


Re, 


Re 


k  1 

9  \ 9 


Units 


215  for  location 
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INPUT  TABLE  1.  Continued 


Card 

No.  Columns  Format 


Units 


kJ 

(Conti nued) 


21-25 


26-30 


31-35 


36  -40 


|  0.1B‘  +  C(0.9  +  0.11CARB)(1  +  0.25B*  )(pe/pW^ 

)  LORN  cone  criterion 

t  ) 

■g  V  7  —  Fully  turbulent 

Positive  NREYCR  -  Transitional  heating 

Negative  NREYCR  -  No  transitional 
heating  (abrupt  transition) 

IRON  —  Body  angle  definition  flag 

-1  —  Circle  curve  fit 

0  —  Circle  fit  for  laminar  flow,  angle 
averaging  for  turbulent  flow  with 
damping  logic  for  negative  curvature 
(recommended  for  shape  change) 

1  —  Angle  averaging 

1CARB  —  Carbon  flag  in  PANT  transition 
cri ter ion  (NREYCR  =  &) 

0  —  Uses  0.9  factor  in  *  calculation 
(recommended  for  carbon  materials! 

1  —  Uses  1.0  factor  in  y  calculation 
1ATH  —  Atmosphere- type  flag  for  LG  -  1 

1  --  Standard  kwajalein  atmosphere 

2  —  U.S.  standard  atmosphere,  1962 

.  1ATMS  —  Atmosphere- type  flag  for  LG  -  -1 


1  --  No  atmosphere 

2  —  1962  standard 

3  —  15  deg  north  annual 

4  --  30  deg  ncrth  January 
6  —  30  deg  ncrth  July 

6  --  45  deg  north  January 

7  —  45  deg  north  July 

8  —  45  deg  north  spring/fall 

9  —  6U  deg  north  January 

10  —  60  deg  north  January  (cold) 
U  —  60  deg  north  January  (warm) 
\2  --  60  deg  north  July 
\3  --  75  deg  north  January 
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*  f 
f* 


INPUT  TABLE  1.  Continued 


Card 

No.  Columns  Format  Data  Units 


14  --  75  dey  north  January  (cold) 

15  —  75  deg  north  January  (warm) 

16  —  76  deg  north  July 

17  —  Tonapah  winter 

18  --  Tonapah  spring 

19  --  Tonapah  summer 
2U  --  Tonapah  fall 


*3  41-45  15  I MOD  —  Nose  shape  modification.  The  code 

(Continued)  does  not  allow  the  surface  angle  between 

two  adjacent  points  to  exceed  89.5°.  The 
user  has  the  option  to  specify  whether  the 
point  (I  +  1)  or  (I  -  1)  is  used  to  control 
the  recession  at  point  I.  IMOD  -  1  usually 
results  in  larger  tiinesteps  in  severe 
erosion  environments. 


0  —  Use  I  and  (I  +  1) 
1  --  Use  I  and  (I  -  1) 


1-  2 

12 

Nonzero  entry  for  the  last  card  in  the 
table 

3-14 

E12.5 

01 TT  --  Environment  time  step  (at  the  end 
of  each  interval  new  environment  will  be 

calculated) 

Sec 

15-26 

E12.5 

TIM  --  Time  for  change  of  DLTT 

Sec 

Same  as 

Card  4  for 

various  output  intervals 

(etc.) 


(Initial  tlfli*) 

Tl  TIM 


TIM 


T]M 


(final  time) 
Ti 


|— X— K-K-f 


OLT  T 


Trajectory  calculation  option  only 
(Read  only  if  L6  =  -1) 
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INPUT  TABLE  1.  Concluded 


Card 

No. 

Columns 

Format 

Datt 

Units 

*6 

3-14 

E12.5 

HTJI  --  Altitude  at  which  trajectory 
calculations  ;tart 

Ft 

15-26 

E12.5 

HSCI  —  Altitude  at  which  shape  change 
calculations  start  (typically  shape 
change  should  be  limited  to  altitudes 
below  150,000  ft) 

Ft 

27-38 

E12.5 

HSCF  —  Altitude  at  which  shape  change 
circulations  terminate 

Ft 

*7 

39-50 

E12.5 

HTJF  —  Altitude  at  which  trajectory 
calculations  terminate 

Ft 

3-14 

E12.5 

TUT  —  Initial  longitude  for  trajectory 
calculation 

Deg 

15-26 

E12.5 

PHGO  —  Initial  latitude  for  trajectory 
calculation 

Deg 

27-38 

E12.5 

VELTJ  —  Initial  velocity  for  trajectory 
calculation 

Ft/sec 

39-50 

E12.5 

GAHGO  --  Initial  flight  path  angle 
(should  be  negative  value  for  reentry 
trajectory) 

Deg 

51-62 

E12.5 

SIGGD  —  Initial  heading 

Deg 
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INPUT  TABLE  2.  ENVIRONMENT  TABLE 


This  table  inputs  environment  conditions  according  to  L6  flag  of 


I  nput 

Table  1. 

Maximum  of  60  entries  are  allowed  in  thi ; 

table. 

1.  Input  for  Flight  Environment,  LG  : 

Card 

No. 

Columns 

Format 

Oat.i 

Units 

1 

1-  2 

12 

Enter  02  (table  number) 

-- 

2 

1-  2 

12 

Nonzero  entry  for  the  last  card 
table 

in  the 

3-14 

E12.5 

Time 

Sec 

15-26 

E12.5 

Altitude 

Ft 

27-38 

E12.5 

Velocity 

Ft/sec 

3 

Same  as 

Card  2  for 

increasing  time 

(etc. 

) 

2.  Input 

for  Wind  Tunnel  Environment,  LG 

2-1 

Card 

No. 

Columns 

Format 

Data 

Units 

1 

1-  2 

12 

Enter  02  (table  number) 

-- 

2 

1-  2 

12 

Nonzero  entry  for  the  last  card 
table 

in  the 

3-14 

E12.5 

Time 

Sec 

15-26 

E12.5 

Free  stream  total  pressure 

psia 

27-38 

E12.5 

Free  stream  total  temperature 

°F 

39-50 

E12.5 

Free  stream  Mach  No.  (a  constant 
on  Card  2  only) 

. ,  entered 

3 

Same  as 

Card  2  for 

increasing  time 

(etc. ) 
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Environment  down  the  range  Projectile  history 


INPUT  TABLE  2.  Continued 


3.  Input  for  Ballistic  Range  Environment,  LG  c  3 


Card 

No. 

Columns 

Format 

Date 

Units 

1 

1-  2 

12 

Enter  02  (table  number ) 

-- 

2 

1-  2 

12 

Nonzero  entry  for  the  last  card  in 
table 

the 

-- 

3-14 

E12.5 

Time 

Sec 

15-26 

E12.5 

Projectile  displacemert  down  the  range 

Ft 

27-38 

E12.5 

Vfcloci ty 

Ft/ sec 

(3+n) 

Same  as 

Card  2  for 

increasing  time 

( n+1 ) 

1-  2 

12 

Nonzero  entry  for  the  last  card  in 
table 

the 

— 

3-14 

E12.5 

Distance 

Ft 

15-26 

E12.5 

Pressure 

Atm 

27-38 

E12.5 

Temperature  (Card  (n+1)  only) 

°R 

( n+2) 

Same  as 

Card  (n+1) 

with  increasing  distance 

(etc.) 

4.  Input 

for  General  Environment,  LG  *  4 

Card 

No. 

Columns 

Format 

Data 

Units 

1 

1-  2 

12 

Enter  02  (table  number) 

-- 

2 

1-  2 

12 

Nonzero  entry  for  the  last  card  in 
table 

the 

-- 

3-14 

E12.5 

Time 

Sec 

15-26 

E12.5 

Free  stream  static  pressure 

Atm 

27-38 

E12.5 

Free  stream  static  temperature 

°R 

39-50 

E12.5 

Free  stream  velocity 

Ft/sec 

3 

Same  as 

Card  2  for 

Increasing  time 

(etc.) 
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INPUT  TABLE  2.  Concluded 
5.  Input  for  Arc  Heater  Environment,  LG  =  5 


Card 

No. 

Columns 

Format 

Data 

Units 

1 

1-  2 

12 

Enter  02  (table  number) 

-- 

2 

1-  2 

12 

Nonzero  entry  for  the  last  card  in 
displacement  history  subset 

the 

-- 

3-14 

E12.5 

Time 

Sec 

15-26 

E12.5 

Displacement 

Inch 

(3+n) 

Same  as 

Card  2  for 

increasiny  time 

(n+1) 

1-  2 

12 

Nonzero  entry  for  the  last  card  in 
flow  conditions  subset 

the 

-- 

3-14 

£12.5 

Time 

Sec 

15-26 

E12.5 

Free  stream  total  pressure  (pt^) 

Atm 

27-38 

E12.5 

Free  stream  total  enthalpy  (ht^) 

Btu/1 bm 

i+2-*m) 

Same  as 

Card  (n+1) 

with  increasing  time 

(m+1) 

1-  2 

12 

Nonzero  entry  for  the  last  card  in 
table 

the 

-- 

3-14 

E12.5 

Distance  from  nozzle  exit 

Inch 

15-26 

E12.5 

Normal  shock  total  pressure  ratio 

(Ptg/Ptj) 

-- 

(m+2) 

Same  as 

Card  m+1  with  increasiny  distance 
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INPUT  TABLE  3.  INITIAL  CONFIGURATION  AND  IN-DEPTH  CONDUCTION  GRID  PARAMETERS 


Card 


Columns 

Format 

Data 

Units 

1-  2 

12 

Enter  03  (table  number) 

-- 

1-  5 

15 

NS  --  Number  of  points  on  the  heated 
surface  of  the  body  (maximum  TO  points) 

-- 

>0  --  Sphere-core  shape  option 

(applicable  only  for  single 
material  todies) 

<0  —  General  shape  option 

6-10 

15 

NPN  --  Number  of  points  on  the  nose; 
applicable  only  to  sphere-cone  option 
(NS  >  0) 

_  _ 

11-15 

15 

MAT  --  Material  index  for  single  material 
nosetip.  If  the  nosetip  is  multimaterial 
(maximum  of  six  in-depth  materials, 
general  shape  option  only)  it  may  be 
entered  as  zero. 

1-  2 

Blank 

3-14 

E12.5 

RNI  --  Initial  nose  radius* 

Inch 

15-26 

E12.5 

ZMAX  —  Maximum  axial  length  (required 
input  for  sphere-cone  option  only) 

Inch 

ZMAX  =  7.-coordi nate  of  the  last  point  on 
the  sphere-cone 

27-38 

E12.5 

THETA  --  Initial  cone  half  angle  (required 
input  for  sphere-cone  option  only) 

Deg 

39-50 

E12.5 

OX  —  Axial  position  of  the  origin  of  the 
rays 

Inch 

<0  --  Flat  back  option 
>0  --  Plug  option 


^Used  to  estimate  transition  altitude  (or  time),  and  for  scaling  of  body 
input  information 
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INPUT  TABLE  3.  Continued 


Card 

No.  Col umns  format  Data  Units 


3 

( Continued) 


The  following  sketches  illustrate  the 
nosetip  configuration/location  of  the 
origin  of  the  ray's  combination  which 
are  referred  to  as  flat  back  or  plug 
configurations. 


a.  Flat  back  (OX  =  ZMAX) 


b.  Plug  (OX  <  ZMAX) 


51-62  E12.5  TS  —  Initial  body  temperature.  This  input 

wTl  1  be  overridden  if  surface  temperature 
distribution  is  input  via  Table  07.  °R 

63-74  E12.5  STRD  —  (Transient  option  only,  1SS  =  0.) 

Maximum  surface  temperature  rise  desired 

between  time  steps.  If  it  is  less  than 

49°R  or  greater  than  201°R,  it  is  set  to 

75°R .  °R 


General  Shape  Option 


The  generalized  shape/interface  option  can  be  thought  of  as  describing 
the  boundary  lines  for  each  material  in  the  vehicle.  Each  material  interface 
should  be  described  as  a  closed  loop  (i.e.,  one  point  should  be  specified 
twice).  A  simple  sample  will  best  illustrate  the  use  of  the  new  option. 
Consider  a  vehicle  modeled  by  the  following  geometry: 
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aag-.M*-. 


Table  3  input  for  this  geometry  is  shown  in  Figure  3-1  with  10  surface 


points,  models  the  interface  with  22  input  points,  and  has  1  plug  point. 
Although  the  configuration  appears  as  a  flat  back,  the  user  has  specified  a 
plug  option  (Figure  b.  Page  3-14)  so  that  the  origin  of  rays  may  be  placed 
closer  to  the  nosetip  for  better  implicit  layer  definition.  As  such,  the  user 
must  specify  the  flat  back  face  as  an  unheated  surface.  From  this  data,  the 
cod,.  will  calculate  the  coordinates  and  material  flags  of  the  interface 
intersections  with  ray  within  the  implicit  layer,  and  the  material  flag 
indices  for  the  explicit  grid,  NMAT. 

Card 

No .  Col umns  Format  Data  Units 

4  1-20  F10.3.F8.3.I2  ZSP( I) ,  RSP( I) ,  NB1( I)  —  For  I  =  l,  NS. 

Body  pofnt  coord fnates  and  material 

indices  for  the  surface  of  the  vehicle.  Inch 

Ente;  one  surface  point  per  card. 

Enter  as  many  cards  as  there  are  surface 
points  (NS) . 

5  1-5  15  NIF  --  Number  of  points  used 

to  describe  the  interfaces 


INPUT  TABLE  3.  Continued 

Card 

No.  Col  urn  ns  Format  Data  Units 

6  1-20  F10.3,F8.3,I2  ZIS(I),  RISC  I ) .  NBS( I)  -  For  I  =  1,  NIF.  Inch 

Coordinates  and  material  indices  for  the 
interface  locations. 

Enter  one  surface  point  per  card. 

Enter  as  many  cards  as  there  are 
interface  points  (NI?). 

If  the  geometry  consists  of  only  one  material,  tne  interface  boundary  is 
described  first  by  tracing  the  surface  body  points  then  closing  the  loop. 

---------------Plug  option  only  ---------------- 

(Read  only  If  OX  >  0) 

Plug  points  Identify  the  remaining  surface  points  that  are  not  on  the  heated 
surface.  Note  if  a  flat  back  configuration  is  specified  as  a  plug  (Figure  b, 
page  3-14),  that  is  OX  >  0,  then  a  single  plug  point  is  required,  connecting 
the  last  surface  point  to  the  axis,  which  describes  the  unheated  flat  back 
face. 

Card 

No.  Colimns  Format  Data  Units 

5  1-2  12  NC  --  F1  ag  to  read  the  coordinates  of  the 

Body  points 

*  0  --  Keep  reading 

f  0  -■  Stop  reading.  This  indicates  that 
the  card  is  the  last  of  its  kind. 

3-14  E12.5  ZSP  --  Body  point  axial  length,  (z)  Inch 

15-26  E12.5  RSP  --  Body  point  radial  length,  (r)  Inch 
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INPUT  TABLE  3.  Continued 

--------------  In-depth  grid  setup  ------------ 

(Transient  option  only:  ISS  -  0) 

This  In-depth  grid  definition  and  nomenclature  are  shown  In  the  following 
sketch. 


The  grid  size  and  distribution  are  problem  dependent.  In  general, 
there  is  no  rule  as  to  what  the  optimum  value  of  the  grid  size  is;  and  one  has 
to  perform  some  nunerical  experiments  to  arrive  at  the  optimum  value.  The 
acceptable  solution  to  the  problem  Is  the  one  which  does  not  change  when  the 
grid  size  is  further  refined. 

In  order  to  obtain  a  rough  estimate  of  the  thickness  of  the  surface 
layer,  we  use  the  results  of  steady-state  analysis  of  a  semi-infinite  solid 
with  constant  surface  temperature  (or  heat  flux)  and  recession  rate,  s.  It 
can  be  shown  that  the  thermal  penetration  depth  in  t»—  solid  is 

Dp  ■  2 .3q 
s 


46 


INPUT  TABLE  3.  Continued 


where  a  is  the  material  thermal  diffusivity  and  Dp  is  defined  to  be  the 
distance  from  the  receding  surface  to  where  the  temperature  drops  to 
10  percent  of  the  surface  temperature. 

In  the  nosetip  application,  if  we  can  estimate  a  characteristic 
recession  rate,  s,  we  may  state  that  the  surfac*;  layer  thickness  should  be 
greater  than  or  at  least  equal  to  Dp  obtained  f tom  the  above  formula. 

For  plug  configuration  the  position  of  the  origin  of  the  rays  on  the 
axis  of  symmetry  is  input  by  the  user.  As  a  guide  to  determine  the  position 
of  this  origin,  it  should  be  noted  that  for  accuracy  of  computations  we  desire 
the  rays  to  be  as  close  to  the  surface  normals  .is  possible.  On  the  other 
hand,  the  distance  OX  should  be  large  enough  to  allow  the  surface  to  recede 
without  getting  too  close  to  the  origin  of  the  rays.  The  computations  are  set 
to  stop  if  the  distance  from  the  back  of  the  surface  layer  to  the  origin  of 
the  rays  is  anywhere  smaller  than  DSMOVE  unless  the  automatic  shifting  option 
is  specified.  In  the  former  case,  the  computations  can  be  continued  by 
relocating  the  origiq  of  the  rays  and  using  the  code  restart  capability.  The 
input  format  for  the  in-depth  grid  is  as  follows: 

Implicit  Grid 

Card 


No. 

Columns 

Format 

Data 

Units 

6 

1-  5 

12 

LL  —  Total  number  of  implicit  nodes 
fnodlets)  along  each  ray  (maximum  15) 

— 

7 

1-  2 

12 

Blank 

3-74 

6E12.5 

DELN(I)  —  For  I»2,  LL;  normalized 
nodlet  spacing  (normalized  distance 

— 

between  nodlets,  from  surface  inwards 
must  sum  to  unity;  maximum  of  14). 

If  uniform  spacing  Is  desired,  enter 
only  one  spacing,  DELN(2);  i.e., 

1/ ( LL-1 )  *  uniform  spacing. 


4/ 


Not  input  for 
uniform  grid  spacing 


INPUT  TABLE  3.  Concluded 


Card 

No.  Columns  Format  Data 


8 

1-  2 

12 

B1  ank 

3-74 

6E12.5 

DEL( I )  --  For  1=1,  NS;  surface  layer 
tffiTckness  along  each  ny.  If  uniform 
thickness  is  desired,  enter  only  one 
thickness,  DEL(l) 

9 

1-  5 

15 

IL  --  Number  of  explicit  nodes  in  the 
^direction  (maximum  63) 

6-10 

15 

JL  —  Number  of  explicit  nodes  in  the 
^direction  (maximum  25) 

10 

1-  2 

Bl  ank 

3-74 

6E12.5 

XDIF(I)  —  For  I  =  2,  IL  (six  to  a  card), 
X-direction  distance  between  grid  nodes. 
For  unifom  grid  spacing  in  both  X  and  Y 
directions  enter  one  value  only,  XDIF(2) 

11 

l-  2 

Blank 

3-74 

6E12.5 

YDIF(J)  —  For  J  =  2,  JL  (six  to  a  card), 
Y-direction  distance  between  grid  nodes 

Card 

-  -  -  Trajectory  calculation  option  only  -  -  -  -  -  - 

(Read  only  if  LG  -  -1) 

No. 

Col umns 

Format 

Data 

*12 

3-14 

E12.5 

THETAC  —  Frustum  angle  to  be  used  in 
calculating  the  aft  body  drag 

15-26 

E12.5 

RBAS  --  Vehicle  base  radius 

27-38 

E12.5 

VLN  —  Vehicle  axial  length 

39-50 

E12.5 

XCG  —  Vehicle  center  of  gravity  location 
from  stagnation  point 

51-62 

E12.5 

WT  —  Vehicle  weight 

Units 


Inch 


Inch 

Inch 


Units 

Deg 

Inch 

Inch 

Inch 

lbm 
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INPUT  TABLE  4.  TRANSITION  TABLE 


Enter  this  table  only  when  NREYCR  *  2,  3,  or  4.  Maximum  of  30  entries 
are  allowed  in  this  table. 


1.  Input  tor  Re@  versus  Me,  NREYCR  =  2 


Card 

No. 

Columns 

Format 

Data 

Units 

I 

1-  2 

12 

Enter  04  (table  nunber) 

-- 

2 

1-  2 

12 

Nonzero  entry  for  the  last  card  in  the 
table 

-- 

3-14 

E12.5 

AM  —  Local  edge  Mach  No.  (Me) 

-- 

15-26 

E12.5 

REM  —  Critical  momentum  thickness 

Reynolds  number  (Re0) 

-- 

3 

Same  as 

Card  2  for 

increasing  Me 

(etc.) 

2.  Input  for  Res  versus  Mg,  NREYCR  =  3 

1 

1-  2 

12 

Enter  04  (table  number) 

-- 

2 

1-  2 

12 

Nonzero  entry  for  the  last  card  in  the 
table 

-- 

3-14 

E12.5 

AM  —  Local  edge  Mach  No.  (Mg) 

-- 

15-26 

E12.5 * 

REM  —  Critical  steam-length  Reynolds 
number  (Res) 

-- 

3 

Same  as 

Card  2  for 

Increasing  Mg 

(etc.) 
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INPUT  TABLE  4.  Concluded 


Columns  Format 


Data 


Units 


3.  Input  for  axial  transition  location 
versus  altitude  (time),  I  IRE  VCR  =  4 


1-  2 

12 

Enter  04  (table  number) 

1-  2 

12 

Nonzero  entry  for  the  last  carU 
table 

in  the 

3-14 

E12.5 

AM  —  Altitude  (LG  =  1  or  3)  or 
ILG  ■  2,  4,  or  5) 

time 

15-26 

E12.5 

REM  --  Axial  transition  location 
current  stagnation  point 

from 

Same  as 

Card  2  for 

increasing  altitude  (time) 

Ft  or 
sec 


Inch 
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INPUT  TA8LE  5.  PARTICLE  ENVIRONMENT 


This  table  is  not  required  for  clear  air  calculations.  Three  options 
are  allowed  for  cloud  characterization  with  a  maximum  of  50  entries  in  each. 


Card 

No. 

1 

*2 


Columns  Format 


1-  2 

1-  5 


6-10 


11-15 


16-20 


1-  2 
3-10 


12 

15 


15 


15 


15 


12 

F8.1 


Date 

Enter  05  (table  nunbei ) 

NOSLO  --  Shock  layer  particle  interaction 
flag- 

<0  _  Interaction  model  is  a  function  of 

the  altitude  and  is  input  in  the 
following  tables 

1  --  No  interaction 

2  —  Jaffe  model 

3  _  Reinecke  and  Waldman  model 

4  —  Ranger  and  Nicholls  model 
NOHEAL  —  Ablation  healing  flag 

0  --  No  crater  healing 

1  —  Crater  healing 

METER  —  Input  units  flag 

0  —  Altitudes  in  feet 

1  --  Altitudes  in  meters 

INPUT  —  Distribution  flag 

0  --  Constant  particle  size 

1  --  Built-in  distributions 

Z  —  Input  distributions 

Option  2  (INPUT  =  0) 

NC  —  Nonzero  entry  for  last  card 

ALCL  —  Altitude  (LG  »  1)  or  distance 
along  range  (LG  1  3) 


Units 


Ft  or 
meter 


51 


INPUT  TABLE  5.  Continued 


Card 

No.  Columns  Format 

11-20  F10. 2 

21-30  F10. 2 

31-40  F10.2 

*  41-50  F10. 2 


Oa  ,a 

021  —  Liquid  water  :ontent 

OPI  —  Particle  diamiter 

361  —  Particle  specific  gravity 

TYP  --  Particle  type 

1.0  —  Ice 

2.0  —  Small  snow 

3.0  —  Large  snow 


4.0  --  Rain 

60  II  NDEM  --  Particle/shock-layer  interaction 

model.  Not  used  if  NOSLO  >  0. 

1  --  No  interaction 

2  —  Jaffe  model 

3  —  Reinecke  and  Waldman  model 

4  —  Ranger  and  Nicholls  model 

4  Same  as  Card  3  for  increasing  or  decreasing  altitude 


Units 
gm/ni^ 
Mi  crons 


(etc.) 

Option  2  (INPUT  =  1) 

*3  1-  2 

12 

NC 

I 

3-10 

F8.1 

ALC1 

11-20 

F10.2 

CZJ. 

• 

Same  as  Option  1 

21-30 

F10.2 

TYP  j 

r  «■» 

31-40 

F10.2 

SGJ.  ] 

i 

50 

11 

NOEM 

) 
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INPUT  TABLE  5.  Concluded 


f  Card 

t'i  fcl  _ 


Columns 

Format 

Data 

Uni  ts 

Option  3  (INPUT  =  2) 

i-  2 

12 

NC 

- 

3-10 

F8.1 

ALCL 

( 

_ 

i  Same  as  ;.>pt  on  1 

11-20 

F10.2 

CZI 

- 

21-30 

F10.2 

TYP 

/ 

- 

31-40 

F10.2 

DMAX 

--  Maximum  particular  diameter 

Mi  crons 

41-50 

F10.2 

ANO  - 

•-  'i,.  I  Parameters  used  to  describe 

(  distribution  function 

?  -Xdp 

51-60 

F10.2 

ALAM 

-  1  N(D)  =  N^e 

- 

61-70 

F10.2 

SGI 

}  Same  as  Option  1 

- 

80 

12 

NDEM 

) 
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INPUT  TABLE  6.  MATERIAL  PROPERTIES 


Maximum  of  three  material  inputs  are  allowed  in  this  table.  Each 
material  input  is  treated  according  to  the  following  forma;.. 


Card 

No. 


Columns  Format 


Data 


11-15 


1-  2 

3-14 

15-26 

27-36 


15 


E12.5 

E12.5 

E12.5 


0  —  No  erosion 

1  —  Generalized  G-law  erosion 

2  --  Carbon-phenolic  law 

3  —  Tungsten  law 

JROUGH  --  Laminar  heating  augmentation  flag 
0  —  No  augmentation 

1  —  Transition  proximity  augmentation 

2  --  Transition  proximity  and  particle 

stirring  augmentation  with  correla¬ 
tion  of  metal  lie  data 

H  .  0.1193  p  u  ( .aJBiast  J>  ) 

3  --  Particle  stirring  augmentation  with 

correlation  for  graphite  data 


0.098 


» . Impact  c 


(1  ♦  G) 


Blank 

RUFL  --  Intrinsic  roughness  height 

RUFMAX  —  Maximum  turbulent  roughness 
height  (kt  ) 
max 

RUF1  --  Constant  kj 

>0  —  White-Grabow  scallop  law  is  used: 

=  RUFMAX 


Units 


1 

1-  2 

12 

Enter 

06  (table  nunibe1 

') 

- 

2 

1-  5 

15 

MAT  — 

Material  index 

fer  following  table 

- 

6-10 

15 

NERGDE 

--  Erosion  law 

material  flag 

- 

Mil 

Mil 

MIL-psi0 


kf  =  kjPe"^*^  with  kt 


max 
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INPUT  TABLE  6.  Continued 


Card 


. "l umns 

Format 

Ha. a 

Units 

-0  - 

-  Constant  jrbu  ent  roughness  of 

RUFMAX  is  used 

1-  2 

Blank.  If  NEROUE  =  <i,  this  card  is 

r.ct 

required. 

3-14 

E12.5 

*  > 

Constants  for  generalized  G-law 
i  erosion  (this  card  required  only 

15-26 

E12.5 

B  1 

1  if  NERQDE  =  1)  where: 

-- 

27-38 

E12.5 

c 

G  =  A(up)B  (dp)C  (mp)D  (sin  e)E  (Tw)f 

-- 

>  Up  =  Impact  velocity  in  ft/sec 

-- 

39-50 

E12.5 

0 

dp  -  Particle  diameter  in  micron 

... 

51-62 

E12.5 

£ 

1  mp  =  Mass  of  the  particle  in  lbm 

-- 

63-74 

E12.5 

1 

F  ■ 

Tw  =  Surface  temperature  in  °R 

— 

75-80 

E6.0 

OPAR l H 

>0  - 

-  Minimum-particle  diameter  used  in 

Mi  crons 

carbon-carbon  G-law 

<0  - 

-  Use  graphite  G-law 

Parameters  in 

Blowing  Correction  to  Transfer  Coefficients 

1-  2 

12 

NC  - 

-  Flag,  zero  for  ISS  =  0  or  1.  For 

T?S 

=  2,  use  -1  for  terminal  card  of  an 

intermediate  material  table  and  +1  for 

1  ast 

material  table.  The  input  of  the 

rest 

of  Input  Table  6,  except  card  9  is 

not 

required  if  ISS  =  2. 

-- 

Default 

Values 

3-14 

E12.5 

BLS 

--  Laminar  shear  parameter  0.5 

-- 

15-26 

E12.5 

BLH 

--  Laminar  heating  parameter  0.5 

27-38 

E12.5 

BTS 

--  Turbulent  shear  parameter  0.35 

-- 

39-50 

E12.5 

BTH 

--  Turbulent  heating 

parameter  0.35 

-- 

If  zero  is  entered,  default  values  are 
used. 
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INPUT  TABLE  6.  Continued 


Card 


No. 

Columns 

Format 

Data 

Units 

6 

1-  2 

Blank 

-- 

3-14 

E12.5 

RHO  --  Material  density 

lbm/f 

15-26 

E12.5 

TFO  --  Datum  temperature  for  heat  of 
formati on.  For  JANAF  data  TFO  =  36°R. 

°R 

27-38 

E12.5 

HFO  —  Heat  of  formation 

Btu/1  bin 

7 

1-  2 

Blank;  if  non-melting  material,  enter 
a  blank  card 

3-'4 

E12.5 

XLATHT  --  Latent  heat  of  fusion  for 
melting  material 

Btu/lbm 

15-76 

E12.5 

TMELT  --  Melt  temperature 

°R 

27-38 

E12.5 

DTMELT  --  Temperature  difference  over 
which  melt  is  to  occur 

°R 

39-50 

E12.5 

CPSOLD  --  Specific  heat  of  so' 

Btu/1 bm-°R 

51-62 

E12.5 

CPLIQD  -  Specific  heat  of  '  a< 

iwur 

Btu/1bm-°R 

8 

1-  2 

12 

NC  --  Flag,  nominally  zero,  +1  i. 

Terminal  card  of  last  material  pr,  -ty 
table,  -1  marks  terminal  card  of 
intermediate  material  property  tables 

-- 

3-14 

E12.5 

Temperature  (independent  variable) 

°R 

15-26 

E12.5 

Specific  heat 

6tu/1bm-°R 

27-38 

E12.5 

Thermal  conductivity  (not  required  for 
the  steady-state  option) 

Btu/ft-sec-°R 

39-60 

E12.5 

Emissi vity 

... 

Additional  materials  (if  any)  follow  in  a  similar  manner, 
starting  with  Card  2. 
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INPUT  TABLE  6.  Concluded 

Card 


No. 

Columns 

Format 

Data 

Contact  Resistarce 

Units 

9 

1-  5 

15 

f 'ATI  --  Material  numter 

-- 

6-10 

15 

MAT2  --  Material  numter 

-- 

11-22 

E12.5 

VALUE  --  Contact  resistance  between  MAT1 
and  MAT2 

Ft2-S-°R/Btu 

Enter 

contact 

resistance 

between  each  pair  of  materials  only  once,  i 

. e .  ,  if 

MAT1  =  2  and  MAT2  =  5  it  Is  not  necessary  to  enter  a  contact  resistance  for  the 
pair  MAT1  =  5  and  MAT2  =  2.  Terminate  entries  with  a  blank  card;  if  there  are  no 
contact  resistances,  the  blank  card  is  still  required. 
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INPUT  TABLE  7.  SURFACE  OATA 


The  BRLASCC  code  includes  sophisticated  techniques  to  predict  surface 
pressure,  temperature,  and  blowing  rate  di stri jutions  of  ablating  surfaces. 

The  user,  however,  nay  bypass  these  computations  and  enter  the  desired  surface 


pressure,  temperature,  or  blowing  rate  distributions  via  this  table.  Entries 
for  each  variable  are  described  in  the  following  subsections.  Any  combination 
of  one  or  three  subsections  can  be  input.  Maximum  of  60  entries  are  allowed 
in  each  subtable. 


Card 

No.  Columns  Format  Data 


1  1-  2 

2  1-  5 
6-10 


Pressure 

12  Enter  07  (table  number) 

15  1TABL  --  Enter  1  (pressure  table) 

15  IPRFLG 


1  --  Surface  pressure  ratio  given  as 
function  of  axial  distance 


3 


4 


1-  2 

3-14 


15-26 


2  --  Surface  pressure  ratio  given  as 
function  of  dimensionless  stream 
length  (S/R^)  where  Rn  is  the 
nose  radius 


12  Nonzero  entry  for  the  last  card  in  the 

table  (-1,  if  other  surface  tables  follow, 
+1  for  last  surface  table) 

E12.5  Axial  distance  from  staoration  point 

(IPRFLG  =  1) 

or,  dimensionless  stream  length  (S/Rm) 
(IPRFLG  =2) 

E12.5  Pressure  ratio  (p/p^) 


Same  as  Card  3  for  increasing  distance  or  S/R^ 


Units 


Inch 
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INPUT  TABLE  7.  Concluded 


Card 

No. 

m 


Temperature 


m+1 


(m+2) 

(etc.) 


n+1 


Columns  Format 


Data 


1-  5 
6-1015 


1-  2 


3-14 


15-26 
Same  as 

1-  5 
6-10 


15  I T ABL  --  Enter  2  (temperature  table) 

I  THE LG 

1  --  Surface  temperature  as  a  function  of 

axial  distance 

2  --  Surface  temperature  as  a  function  of 

dimensionless  scream  length  (S/R^) 

12  Nonzero  entry  for  last  card  in  table 

(-1,  if  other  surface  tables  follow, 

+1  for  last  surface  table) 

E12.5  Axial  distance  from  stagnation  point 

( ITMFLG  =  1) 

or  dimensionless  stream  length  (S/Rjj) 
(ITMFLG  =2) 

E12.5  Temperature 

Card  m+1  for  increasing  distance 

Blowing  Rate 

ITABL  --  Enter  3  (blowing  rate  table) 


15 

15 


(n+2) 

(etc.) 


1-  2 

3-14 

15-26 
Same  as 


IBRFLG 

1  --  Blowing  rate  (B*)  as  a  function  of 

axial  distance 

2  —  Blowing  rate  (B‘)  as  a  function  of 

dimensionless  stream  length  (S/Rf)) 

Nonzero  entry  for  last  card  in  table 
(-1,  if  other  surface  tables  follow, 

+1  for  last  surface  table) 

Axial  distance  from  stagnation  point 
(IBRFLG  =  1) 

or  dimensionless  stream  length  (S/Rfj). 
(IBRFLG  =  2) 


E12.5  Blowing  parameter  B'  h  Pwvw/peue^M 
Card  n+1  for  increasing  distance 


12 


E12.5 


Units 


Inch 


°R 


Inch 
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INPUT  TABLE  8.  SHOCK  SHAPE  DATA 


The  ASC  code  includes  state-of-the-art  shock  shape  nredici.lon 
capability.  The  user  may,  however,  bypass  the  code  shocK  ‘>ape  calculations 
and  input  the  shock  coordinates  and  shock  angle  .is  a  function  of  shock  radial 
coordinate  via  this  table.  Maximum  r  150  entries  are  allowed  in  this  table. 
Card 


No. 

Columns 

Format 

Data 

Units 

1 

1-  2 

12 

Enter  08  (tab. a  number i 

— 

2 

1-  5 

15 

ISHFLG  --  Shock  shape  flag 

-- 

1  —  Shock  angle  given  as  function  of 
y-coordinate 

2  —  Shock  angle  given  as  function  of 
dimensionless  y-coordinate  (y/Rfl), 
where  Rfl  is  the  nose  radius 

3 

1-  2 

12 

Nonzero  entry  for  the  last  card  in  the 
table 

-- 

3-14 

E12.5 

y-coordinate  (for  ISHFLG  =  1)  or 
dimensionless  y-coordinate  (y/Rn), 

(ISHFLG  =  2) 

Inch 

15-26 

E12.5 

Shock  angle 

Oegrees 

If  Table  5  is  entered  with  the  parameter  NOSLO  *  2,  3,  or  4, 
the  x-coordinate1"  of  the  shock  must  also  be  entered 
according  to  the  following  format: 

27-38 

E12.5 

x-coordinat°  (for  ISHFLG  =  1)  or 
dimensionless  x-coordinate  (x/Rr) 

(ISHFLG  =  2) 

Inch 

4 

Same  as 

Card  2  for 

Increasing  y  or  y/Rfl 

For  normal  shock  calculation,  make  a  table  having  one  entry 
with  shock  angle  of  90°. 


^Thi s  input  is  required  for  weather  calculations  with  demise. 


INPUT  TABLE  9.  SURFACE  THERMOCHEMISTRY 


Table  9  consists  of  the  parameters  necessary  for  the  surface  energy 
balance  formulation.  For  multimateri al  nosetips,  the  same  format  is  repeated 
for  each  material.  Different  material  thermcch jmi stry  tables  must  be 
separated  by  a  blank  card,  and  the  final  table  is  terminated  by  two  blank 


cards. 

Card 

No. 

Columns 

Format 

Data 

Units 

1 

1-  2 

12 

Enter  09  (table 

number) 

-- 

2 

1-  5 

15 

MAT  --  Material 

index  number 

-- 

3 

1-  2 

Blank 

-- 

3-14 

E12.5 

CMH  —  Ratio  of 

mass  to  heat  transfer 

coefficients  (typically  1.0) 


Edge  Tables 

If  diffusion  coefficients  are  not  equal  or  if  the  ratio  Cm/Ch  is  not 
unity,  then  the  surface  energy  balance  requires  data  about  the  edge  gases  of 
the  boundary  layer.  These  data  are  provided  in  special  "edge  tables"  which 
precede  each  pressure  section  of  the  surface  tables  (the  various  sections  of 
the  surface  tables  are  described  below).  The  independent  variables  for  an 
edge  table  are  pressure  and  temperature.  Dependent  variables  are  hpw  and 


*  T 


the 

sum  Ez-jeh-|w. 

Card 

No. 

Columns 

Format 

Data 

Units 

4 

1-  8 

F8.5 

Pressure 

Atm 

9-24 

Blank 

1 

-- 

25-33 

F9.4 

Temperature 

(°R  if  negative 
in  which  case 
enthalpies  below 
are  Btu/lbm) 
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INPUT  TABLE  9.  Continued 


Card 

No. 

Columns 

Format 

Data 

Units 

4 

34-38 

F5.3 

Unequal  diffusion  exponent  y 

(Continued) 

*  T 

39-47 

F9.3 

Summation  £ziehiw 

cal/gr 
(Btu/lbm  if 
temperature 
is  entered 
with  minus 
sign) 

T 

48-56 

F9.3 

Enthalpy  of  edge  gases  hew 
evaluated  at  Tw 

cal/gr 
(Btu/lbm  if 
temperature 
is  entered 
with  minus 
sign) 

57-58 

12 

-1  (flag  signifying  that  this  card 
is  part  of  the  edge  gas  table) 

••  «• 

5  Same  as  Card  4  for  remaining  entries  in  "edge  table"  for  this 
pressure,  maximum  of  12  temperatures  for  each  pressure. 

(etc. ) 


The  table  length  is  limited  to  5  pressure  sets  (it  may  have  only 
1  pressure  set)  with  not  more  than  12  nor  less  than  3  temperature  entries  in 
each  set.  The  series  of  temperature  values  may  be  different  for  the  edge 
table  at  each  pressure  set.  The  table  Is  organized  as  a  series  of  sections, 
each  representing  one  pressure  and  each  preceding  the  corresponding  pressure 
group  of  the  surface  thermochemistry  deck  as  described  below.  The  temperature 
entries  within  each  section  must  be  ordered,  either  ascending  or  descending. 
Similarly,  the  pressures  must  be  ordered  either  ascending  or  descending. 

Decks  generated  by  Aerotherm  thermochemistry  programs  will  have  been 
automatically  ordered  properly. 
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INPUT  TABLE  9.  Continued 


Surface  Thermochemistry  Tables 
Description  of  Surface  Thermochemical  Tables 

This  table  comprises  a  series  of  sections.  Each  section  represents  one 

pressure  and  one  transfer  coefficient  value.  More  than  one  transfer 

coefficient  may  be  necessary  if  the  effects  of  kinetics  on  the  surface 

response  are  considered.  Nondimensional  ablation  rate,  B'  ,  forms  the  third 

tc 

independent  variable  within  a  given  section.  The  table  has  three  dependent 
variables:  £ziwh^w,  hww,  and  Tw. 

The  Aerotherm  thermochemistry  codes  generate  separate  groups  for  each 
pressure,  one  at  a  time.  All  these  groups  together  make  up  the  surface 
thermochemistry  deck.  Within  each  pressure  group  the  transfer  coefficient 
values  will  be  ordered.  Within  each  transfer  coefficient  section,  ablation 
rate  entries  need  not  be  ordered  in  any  particular  way  on  the  ablation  rates; 
any  necessary  ordering  is  made  automatically  by  the  cede  as  it  reads  in  the 
data . 

Users  providing  their  own  thermochemistry  decks  must  ensure  that  the 
transfer  coefficients  are  ordered,  but  the  ordering  may  be  either  ascending  or 
descending  in  each  case.  The  surface  thermochemistry  cards  are  identified  by 
a  unity  flat  in  Column  58,  as  described  in  the  format  specification  below. 

The  number  of  pressure  groups  may  not  exceed  5  (and  may  be  only  1);  the 
number  of  transfer  coefficient  values  in  each  pressure  group  may  not  exceed  5 
but  may  be  only  1.  If  no  kinetics  effects  are  to  be  considered,  a  transfer 
coefficient  of  zero  is  acceptable.  The  sequence  of  transfer  coefficient 
values  need  not  be  the  same  in  the  different  pressure  sections.  Within  each 
transfer  coefficient  section  the  number  of  ablation  rate  entries  may  not 
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INPUT  TA8LE  9.  Continued 


exceed  25  and  may  not  be  less  than  2.  The  series  of  ablation  values,  8‘  ,  may 

tc 

be  unique  for  each  section. 

The  °R-Btu/lb  option  described  for  the  edge  cables  may  be  used  for  these 
tables  also. 

Card  Formats 


Columns 

Format 

Data 

Units 

1-  8 

F8.5 

Pressure 

Atm 

9-16 

F8.5 

Transfer  coefficient* 

1 b/ft2-sec 

17-24 

F8.5 

Nondimensional  ablation  rate 
iVpeue^M  =  B* 

-- 

25-33 

F9.4 

Surface  temperature 

°K 

( °R  is  negative 
in  which  case 
enthalpies 
below  are 
Btu/lbn 

34-38 

F5.3 

Unequal  diffusion  exponent  y 

*  T 

39-47 

F9.3 

Summation  Ez i whiw 

ca i /gr 
(Btu/lbm  if 
temperature  Is 
entered  with 
minus  sign) 

48-56 

F9.3 

T 

Enthalpy  of  wall  gases  hww 

cal /gr 
(Btu/lbm  if 
temperature  is 
entered  with 
minus  sign) 

57-58 

12 

Flag  indicating  surface 
tnermochemistry  table  entry 

— 

59-60 

Blank 

_  „ 

*Must  be  provided  by  the-  user. 
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INPUT  TABLE  9.  Continued 


Card 


No. 

Columns 

Format 

Data 

Units 

n  61-66 

( Continued) 

A6 

Chemical  symbol  of  surface  species. 
(Aerothenn  ACE  and  GASKET  prog,  ains  print 
such  symbols  arranged  alphabetically  and 
truncated  from  riyht  end  if  necessary  ' 

67-78 

E12.3 

Nondiinensional  mechanical  fail  rate  = 
m  /peueCM  -  B*  ,  (TBEP) 
fail 

-- 

n+1  Same  as  Card  n  for  remaining  entries  in  th.s  section;  maximum 
of  26  entries  in  each  section. 

Assembled  Thermochemical  Deck 

The  following  sketch  shows  a  picture  of  an  assembled  thermochemical 
data  deck  for  a  given  material  and  severa1  pressures.  The  deck  corresponds 
to  repeating  the  input  described  previously  for  each  pressure. 


s» 


for 

■  Is 

Ow  f  required) 


Sketch  of  surface  thermochemistry  table  makeup  for  a  yiven 
material  Including  leading  constant  cards. 

G5 


INPUT  TABLE  9.  Concluded 


The  surface  equilibrium  data  deck  for  each  material  must  be 
by  a  single  blank  card.  Output  decks  of  Aerotherm  thermochemistry 
not  have  such  a  card,  and  the  user  must  supply  it.  Therriochemi cal 
for  additional  materials  (if  required)  follow  ir  a  similar  manner, 
material  thermochemistry  table  is  terminated  by  two  blank  cards. 
End  of  Input 

The  end  of  the  input  data  deck  for  each  case  is  signaled  by 
card  with  a  -1  punched  in  Columns  1  and  2.  The  end  of  the  overall 
deck  is  signaled  by  a  -1  punched  in  Columns  4  and  5. 


terminated 
programs  do 
data  decks 
The  last 


a  single 
input  data 
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3.2  SAMPLE  PROBLEM 

This  subsection  presents  a  sample  problem  that  demonstrates  the 
capability  of  the  improved  ASC  code,  BRLASCC.  The  problem  consists  of  a  long 
slender  multimaterial  configuration  with  an  in-depth  me’ting  material.  The 
objective  of  the  analysis  is  to  determine  if  tie  outer  material  "burns 
through"  to  the  inner  materials.  Sample  JCL  for  file  assignments  on  VAX/VMS 
systems  and  a  complete  listing  of  the  input  data  is  presented,  while  only  a 
selected  listing  of  the  output  (which  is  voluminous)  is  presented. 

The  vehicle  is  a  slightly  modified  Army  projectile.  The  nose, 
originally  sharp,  has  been  given  a  small  radius  to  facilitate  gridding  in  the 
nose  region.  The  vehicle  is  biconic,  12.5°  on  the  nose  with  an  8°  afterbody. 

The  BRLASCC  solution  failed  to  converge  to  a  proper  surface  energy 
balance  during  Environment  8.  A  restart  was  undertaken  from  Environment  7 
and  is  shown  as  part  of  the  sample  output.  Asterisks  next  to  the  total 
recession  rate  and  B*  tnermochemical  printout  at  Surface  Point  1  on  the  "Body 
Point  Location  and  Surface  Energy  Balance  Results"  summary  page  indicate 
that  nose  modeling  logic  which  prevents  indented  shapes  from  occurring  is 
being  used, 

BRLASCC  indicates  that  burnthrough  does  indeed  occur  at  sometime  near 

0.59  s. 
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The  sample  VAX/VMS  DCL  command  procedure  to  assign  files  for  I/O  and 


run  BRLASCC  is  shown  below. 


$! - ! 

$!  ! 

S!  RUNBRL.COM  ! 

S!  ! 

S!  COMMAND  PROCEDURE  TO  RUN  BRLASCC  ! 

S!  ! 

S! - ! 

S! 

S  SET  NOVERIFY 

$  SET  DEFAULT  [directory  name] 

S! 

S  ASSIGN  RUN.LOG  SYSSPRINT 

S  ASSIGN  RUN.LOG  SYSSOUTPUT 

S! 

$  ASSIGN  input  data  filename  FOR005 
S  ASSIGN  output  data  filename  FOR006 
S! 

$  ASSIGN  ENV1RON.OUT  FOROOl 

S  ASSIGN  RESTART.F14  FOR014 

S  ASSIGN  TRAJECT.DAT  FOR020 

|  RUN  ""“W  «»e]BRLASCC.EXE 

$  ERASE: 

S  RUN _ LOG  -  FSSEARCH(”RUN.LOG”) 

S  IF  RUN _ LOG  EQS-  THEN  GOTO  NO_LOG 

$  DELETE  RUN.LOG;* 

S! 

S  NO _ LOG: 

S  TRAJ  -  FSSEARCHCTRAJECT.DAT”) 

S  IF  TRAJ  .EQS.  THEN  SEXIT 

S  DELETE  TRAJECT.DAT:* 

S! 

S  EXIT 
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SLOPE  BEFORE  POINT  5  HAS  CHANGED  Br  A  FACTOR  OF  20  -  COU’UTE  NEW  ENVIRONMENT 
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make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc . ) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
th  s  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 


Name 

Telephone  Number 
Organization  Address 


FOLD  HERE 


Director 

US  Army  ballistic  Research  Laboratory 
ATTN:  DRSMC-BLA-S  (A) 

Aberdeen  Proving  Ground,  MD  21005 
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